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INTRODUCTION 


The Satellite Telecommunications Analysis and Modeling Program (STAMP), developed 
by Convair, under contract to NASA/Ames provides the user with a flexible and com- 
prehensive tool for the analysis of ITS system requirements. While obtaining minimum 
cost design points, the program enables the user to perform studies over a wide range 
of user requirements and parametric demands. The program utilizes a total system 
approach wherein the ground uplink and downlink, the spacecraft, and the launch 
vehicle are simultaneously synthesized. A steepest descent algorithm is employed 
to determine the minimum total system cost design subject to the fixed user require- 
ments and imposed constraints. In the process of converging to the solution, the 
pertinent subsystem tradeoffs are resolved. 

STAMP provides the system designer with the capability to consider a wide variety of 
communication systems. By appropriate specification of the pertinent input variables, 
the program can be made to cycle through the various types of subsystem elements 
and user requirements to be considered. Minimum cost solutions obtained can be 
compared to determine desired configurations. 

STAMP can handle data in any one or combination of four different data types: audio, 
video, facsimile, and digital data. The program design is such that computations 
are performed on a per beam basis eliminating the need for assuming a nominal area 
of coverage and uniform grouncTsy stemTTo r the~case dTmultiple beams. Each beam 
may individually handle any of the four previously mentioned data types. The logic 
also includes a variety of receive and transmit capability options for the ground 
stations . 

The program can be used as a tool in performing required sensitivity analyses. For 
example, it can be employed to determine the effected change in the total system 
cost resulting from a perturbation of any of the design parameters, the sizing con- 
straints, and/or the subsystem parameter relationships, as well as the key user 
requirements such as S/N, frequency, and area of coverage. 

This report documents STAMP. Section 2 provides a technical analysis 
and a description of the principal techniques employed in the program. 
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TECHNICAL DISCUSSION 


2 . 1 PROGRAM DESCRIPTION 

The ITS synthesis program is essentially structured into five major portions: the 
driver program, the spacecraft system submodels, the uplink and downlink ground 
system models, the steepest descent iterative process, and other required models 
such as the signal propagation and sizing models (see Figure 2-1). 

In its present form it can handle a list of up to 78 constrained parameters including 
spacecraft weight, spacecraft volume, spacecraft prime power supply, and trans- 
mitter power out/ channel at the spacecraft and at the ground uplink facilities. The 
independent design parameters are chosen from the set of ground antenna gains 
(diameters) and receiver noise temperatures (noise figures). 


The synthesis program is completely modularized in the sense that the various 
models have been written as separate subroutines and may be modified or inter- 
changed with a minimum amount of effort. The program essentially consists of 35 
fixed purpose subroutines representing the various models required to synthesize 
a communications system. Figure 2-2 illustrates the subroutine breakdown and their 
main functions. Briefly, the ITS Driver program controls overall execution and 
essentially performs all of the needed initialization. Subroutine EQUAT maintains 
control over the computation of the dependent parameter values . This subroutine 
controls computation of the required transmitter powers and the costs, weight, and 
volume terms for the various subsystems in addition to sizing the power subsystem 
in the spacecraft. It computes the cost of the selected launch vehicle and computes 
total system cost. Subroutine SPCRFT contains the actual logic for the sizing and 
costing of the various spacecraft submodels. As illustrated, each of the various 
submodels is represented as a separate subroutine. Also included is a user's diag- 
nostic routine (DIAGN) to provide specific diagnostic messages during execution. 
These do not include the required general purpose subroutines which are listed in 
Appendix B. 


The general logic flow of the program operation is given in Figure 2-3. Beginning 
with an initial design point and fixed user requirements, the program first computes 
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Figure 2-1 


Synthesis program block diagram 
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ITS DRIVER PROGRAM - CONTROLS OVERALL EXECUTION 


INPTX 

AOC 

DIAM 

PARBDS 

ATTEN 

XBD 

EQUAT 

BNDV 

OPTMUM (MAXBD) 

CHANEL 

CNVRGE 

OUTPT1, OUTPT2 


CONTROLS INPUT 
ORBIT/COVERAGE MODEL 
COMPUTES SPCRFT ANT SIZE 
SIZING MODEL 
SIG. PROP/NOISE MODEL 
CHECKS X Q AGAINST BDS 
COMPUTES Y(X) 

CHECKS Y FOR BD VIOLATION 
MIN. COST CONVERGENCE ROUTINE 
COMPUTES CHANNEL REQUIREMENT 
CONVERGENCE TEST 
CONTROLS OUTPUT 


EQUAT — COMPUTES Y 

GRDFAC 

GDIRCT 

DGDEP 

DISCR 

PXTR 

SPCRFT 

LNCHV 

TCOST 


UPLINK/DOWNLINK MODEL FOR CLASS I/II 
DIRECT RCVR GRD DWNLNK MODEL 
ITER SCHEME ON GRD ANT GAIN AND DIAM 
COMPUTES GRD ANTENNA DISCRIMINATION 
COMPUTES REQ'D XMTR POWER/CHANNEL 
SPCRFT MODEL 
LAUNCH VEHICLE 
TOTAL SYSTEM COST 


SPCRFT (SPACECRAFT SUBMODELS) 


ANTS 

ANTENNA 

STNKP 

ST. KEEPING 

ATTCN 

ATTITUDE CONTROL 

STRUC 

STRUCTURE 

MANPRV 

MANNED PROVISIONS 

THRML 

THERMAL CONTROL 

MTPLXS 

MULTIPLEXER 

TLMCMD 

TELEMETRY & COMMAND 

POWER 

POWER 

XMTR 


RCVRS 

RECEIVERS 

(PLACE) 

SPCRFT TRANSMITTER 

MISCELLANEOUS 

FUNCTIONS 


DIAGN 

DIAGNOSTIC ROUTINE 

EXPL 




QUAD 

PARAMETRIC EQ'S. 



PFUNC 



Figure 2-2. ITS synthesis program subroutine structure. 
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Figure 2-3. Program operation flow diagram. 


the required spacecraft antenna size based on the desired frequency and area of 
coverage. It proceeds to cost the ground systems and to compute the required 
power out/channel via the link equations. The satellite system is sized and its cost 
computed. Following this the total system cost is determined. Parameters are 
checked for boundary violations and the required partial derivatives are found. The 
optimum design parameter perturbation is computed, and, if convergence has not 
yet been achieved, the process is repeated with the newly perturbed design point 
until convergence is achieved. If more cases are to be considered, the program 
repeats the entire process. The program was designed, tested, and debugged using 
the CDC 6400 computing facility at Convair Aerospace and eventually modified for 
use on NASA/Ames IBM 360 system. The IBM version may incorporate the overlay 
feature to conserve required core space. This overlay is described in Appendix C. 
Average core requirement on the CDC system is 76000| words including the routines 
required by the operating system. The IBM 360 version with overlay and use of 
double precision requires approximately 160 kilobytes of storage. The double pre- 
cision is necessary because of the shortened word length on the IBM machine, 32 
bits/word as opposed to 60 bits/word, and the sensitivity of parts of the system to 
small perturbations. 


*76000g = 31744 10 
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2.2 MATHEMATICAL DESIGN OPTIMIZATION TECHNIQUE 

The Convair Aerospace ITS system design optimization program uses a steepest 
descent iterative process, which is a modified gradient technique, to determine the 
minimum cost design point for a given system configuration. At each iteration the 
nominal design point is perturbed in such a way as to cause the maximum rate of 
decrease in total system cost to be consistent with the design constraints. 

The philosophy of the technique is that the direction of steepest descent is given as 
the direction of the negative cost gradient. When a physical constraint (e.g. , satel- 
lite weight) is encountered, the new direction is given by the component of the nega- 
tive cost gradient tangent to the constraint surface. Referring to Figure 2-4, the 
lowest point D in the valley represents the minimum cost design point; the cost con- 
tours represent the increase in cost as the system design deviates from the optimum. 
The constraint boundary represents a "fence" or physical constraint which prohibits 
attainment of the specified goal. The process begins at Point A and proceeds iter- 
atively from one contour to the next to Point B, the first encounter with the constraint. 
It then proceeds along the boundary (in the direction of steepest descent) until it 
finds the lowest Point C. This represents the lowest achievable point consistent 
with the imposed constraint. 



Figure 2-4. Steepest-descent iterative process. 
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If Ay. is the amount by which the j boundary is violated and g represents the grad- 
ient of total system cost, then the perturbation in X, the vector of independent de- 
sign parameters , is given by 

AX = - K [ I -0 (0 T £) 1 0 T 1 g - K 0(0 T 0) Ay. 
cL J y j 

where 0 is the gradient of the constraint relationship. The derivation of this equa- 
tion is given in Appendix A . 

If Ky is set equal to unity, the second term represents, to first order terms, the 
amount by which X must be changed to reduce Ayj to zero. If Ky. is less than one, 
the boundary violation is partially reduced. The first term represents the component 
of the cost gradient parallel to the constraint boundary and is orthogonal to the 
second term. For small AXj (so that the linear approximation is valid) , the first 
term reduces the cost and contributes nothing directly to the boundary violation. If 
no boundaries are violated, then Ay. and 0 are zero, so that the equation reduces 
to 

AX = - K c g 

which represents a perturbation in the direction of the negative of the cost gradient. 

Since K c is arbitrary .some control over the step-size must be introduced. The 
principle used in determining the magnitude of the step is that if successive steps 
have been in the same direction it is desirable to increase the step-size in order to 
speed up the convergence process. On the other hand, if successive steps have 
been in opposite directions, it is desirable to decrease the step-size since the opti- 
mum design point has been overshot. Furthermore, it appears that larger step- 
sizes are most desirable in the early iterations to hasten arrival to the vicinity of 
the optimum design point. Once in the vicinity of the optimum, smaller steps would 
prove more practical. To afford faster convergence, the step size S is replaced with 
a vector of dimension N. This enables the optimization scheme to scale each com- 
ponent relative to its ability to converge. The modification of the scale factor is 
based on the following criteria. Let AX^ and AX 2 be the two previous perturbations 
and AX be the present. 

Furthermore, let = AX^ • AX and C 2 = AX 2 • AX where the indicated products 
are vector dot products . Then 

SCALE = .5 x SCALE if C 1 < 0, 

SCALE = 2 x SCALE if C v C 2 > 0, 

SCALE = 32 if SCALE >32, and 

SCALE = .5 x SCALE if i_> 30 and SCALE > 1 where i = iteration index. 
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With the step-size, S-, thus specified, the appropriate value of K c can be deter- 
mined from 

/ 'A '/ ’ I I _ ■ 

H = I-?(0 T 0) 


K = 
c 


/ o 2 T T -1 
I S. - K y • Ay. (0 0 ) Ay. 


T 

H H 


in the event of a boundary violation, 
mined from _ 


K = 
c 



If the solution is unconstrained K c is deter- 


These are the equations used by the ITS convergence routine. It is important to 
keep S. from growing to a point at which a large number of iterations would be re- 
quired to bring it down. Therefore the upper bound is placed on S^. 

The rate of convergence to the minimum cost design point is dependent upon the 
relative sensitivities of the cost to each of the design parameters. Since the algor- 
ithm operates on the most sensitive parameters, convergence to the optimum value 
of the less sensitive parameters can be slow. For this reason it is often desirable 
to rescale the units of the design parameters to equalize these sensitivities. 

An example of how the process actually converges, consider Figure 2-5 which 
illustrates a particular case using this technique. In this example, a volume con- 
straint of 170 cubic feet has been imposed on the volume of the solar array. Since 
the initial design point violates the constraint, several iterations are required to 
return this quantity to the specified boundary. The process then proceeds along the 
boundary until the minimum cost solution is achieved. Also illustrated is the inter- 
action of the ground downlink antenna size and the ground downlink receiver noise 
temperature during the convergence process. The minimum at iteration three 
actually represents the unconstrained minimum cost design point. The ITS synthesis 
program is designed to provide a brief summary output at each iteration, so that 
the user is made aware of such occurrences. 

2.3 PARAMETRIC SUBSYSTEM MODELS 

The ITS subsystem models include the pertinent elements required to obtain a com- 
plete synthesis of the total system. The cost, weight, and volume parametric re- 
lationships and fitted curves may be found in the Appendix section of Volume 3. 

In general, the curves were fit with one of three types of equations. The forms of 
these equations is given as one of the following: 
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Figure 2-5. Illustration of convergence process. 
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g(z)/lO 

(1) y = 10T K U 

p 

where g(z) = a + bz + cz , z =10 logP, and P represents the value of the 
parameter determining y, 

(2) y = a + bx c , and 

(3) y = a + bx + cx^. 

Other relations are indicated by either linear or exponential fits or are given as 
thruput values. 

The curves used during this study are applicable over a wide range of parametric 
values. Where it was impossible to obtain an accurate fit with a single equation, 
the curve was fit in more than one piece. The breakpoint, in most cases, is vari- 
able and is included in the input lists. 

The system synthesis incorporates a loss model for both the uplink and downlink 
portions of the beam, a noise model, and scale factors to account for redundancy 
in the spacecraft. It contains an orbit and area of coverage model which computes 
the required beamwidths and a separate model to compute the channel characteris- 
tics. Also included is a satellite power subsystem model. 

The loss and noise model is discussed in Section 2.6; the area of coverage model in 
Section 2.5, and the communication links in Section 2.9. A discussion of the re- 
dundancy scale factors is in Section 2.8. The cost elements are in Section 2.4. 

The choice of parameters for the dependent and independent design vectors is dis- 
cussed in Section 2.10, while the power subsystem model is found in Section 2.7. 

2.4 COST MODEL 

The ITS cost model summarized in Figure 2-6 illustrates the breakdown per cost 
element. The model computes the acquisition, non-recurring (development), opera- 
tion, and maintenance cost to form the total cost. 

Table 2-1 is a tabulation of the cost model as implemented in the computer program. 
The associated array name as used in the program is given together with a descrip- 
tive title and whether it is computed from an input curve or is input as a thruput 
value. It should be noted that each of the individual cost items is dimensioned by the 
number of beams. The cost-estimating relationships are given in the Appendix 
section of Volume 3. 

In referencing the table, note that the following definitions are provided: 
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Table 2-1. Cost Model 


I. GROUND UPLINK/DOWNLINK - CLASS I (CLASS II) 


A. Acquisition Cost 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 

11 . 

12 . 

13. 

14. 


Facilities 

Terminal equipment (audio) 

(video) 


ACQC1 (ACQC2) 


Curve 


Transmitter 


(facsimile) 
(digital data) 
(audio) 
(video) 
(facsimile) 
(digital data) 


Antenna 
Receiver 
Standby power 
Test equipment 

Total acquisition cost — 
Class I (Class II) 


TACQC1 (TACQC2) 


1 

Thruput 

ZACQC1 (1,1) 
thru ACQC1 (13,1) 


B. Operations Cost (per year) 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 

11 . 


Terminal equipment (audio) 

(video) 


OPC1 (OPC2) 


Transmitter 


Receiver 

Personnel 


(facsimile) 
(digital data) 
(audio) 
(video) 
(facsimile) 
(digital data) 


Total operation cost — 
Class I (Class n) 


TOPC1 (TOPC2) 


N^oxACQCl (2,1) 
N 2 %xACQC1 (3,1) 
N 3 %xACQC1 (4,1) 
N 4 %xACQC 1 (5,1) 
Curve 


Curve 

N 5 %xACQC 1 (10,1) 
Curve 

N l x ZOPCl (1,1) 
thru OPC1 (10,1) 
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Table 2-1. Cost Model (continued) 
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C. Installation Cost 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 

11 . 

12 . 

13. 


Terminal equipment (audio) 

(video) 


CINST1 (CINST2) 


Transmitter 


(facsimile) 
(digital data) 
(audio) 
(video) 
(facsimile) 
(digital data) 


Antenna 

Receiver 

Standby power 

Design, integration, & Mgmt 

Total installation cost — 
Class I (Class II) 


TINC1 (TINC2) 


Ng%xACQCl (2,1) 
N ? %xACQCl (3,1) 
N g %xACQCl (4,1) 
N 9 %xACQC1 (5,1) 
N 1q %xACQC 1 (6,1) 
N n %xACQCl (7,1) 
N 12 %xACQC1 (8,1) 
N 13 %xACQC 1 (9,1) 
Curve 

N 14 %xACQC1 (11, 1) 
N 15 %xACQC1 (12,1) 
N 16 %x£ACQC 1 (2,1) 
thru ACQC1 (11,1) 

ZCINST1 (1,1) 
thru CINST1 (12,1) 


D. Maintenance Cost (per year) 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 

11 . 

12 . 

13. 


Facility 

Terminal equipment (audio) 

(video) 


CMNT1 (CMNT2) Mj% x ACQC1 (1,1) 


Transmitter 


(facsimile) 
(digital data) 
(audio) 
(video) 
(facsimile) 
(digital data) 


Antenna 
Receiver 
Standby power 

Total maintenance cost — 
Class I (Class II) 


M 2 % x 


M 12 % x 


( 2 , 1 ) 


(12’, I) 


TMNTC1 (TMNTC2) N L x XcMNTl (1,1) 

thru CMNT1 (12,1) 
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Table 2-1. Cost Model (continued) 


II. GROUND DOWNLINK - DIRECT 


A . Unit Recurring 


1 . 

2 . 

3. 

4. 


5. 


6 . 

7. 


Antenna 

Receiver (Beam 1) 
1 (Beam 2) 


4 t * 

Receiver (Beam 6) 


UREC3(1) 

UREC3(2) 

(3) 

(4) 

(5) 

i s (6) 

UREC3(7) 


Curve 


Curve 


B. Acquisition Cost 


1 . 

Antenna 

ACQC3 (1,1) 

UREC3(1) 

2. 

Receiver 

ACQC3 (2,1) 

UREC3(I+1) 

3. 

Total acquisition cost - direct 

TACQC3 

(1) + (2) 

!. Installation Cost 



1 . 

Antenna 

CINST3 (1, 1) 

Curve 

2. 

Receiver 

CINST3 (2,1) 

50% x ACQC3(2,I) 

3. 

Total installation cost - direct 

TACQC3 

(1) + (2) 


D. Maintenance Cost (per year) 

1. Antenna 

2 . Receiver 


CMNT3 (1, 1) 10% x ACQC3(1, 1) 
CMNT3 (2, 1) 10% x ACQC3(2, 1) 

N l x [(1) + (2)] 


3. 


Total maintenance cost - direct 


TMNTC3 



Table 2-1. Cost Model (continued) 
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III. SPACECRAFT SYSTEM 
A. Acquisition Cost 


1 . 

Prime power 

ACQC4(1) 


Curve 

2. 

Secondary power (batteries) 


(2) 




3. 

Power conditioning 






4. 

Power distribution 






5. 

Antenna 






6. 

Transmitter 






7. 

Multiplexer 





f 

8. 

Receiver 




Curve 

9. 

Structure 




Thruput 

10. 

Thermal control 




Curve 

11. 

Stationkeeping 




J 

t 

12. 

Attitude control 




Curve 

13. 

Telemetry and command 




Thruput 

14. 

Manned provisions (life support subsystem) 



Thruput 

15. 

Integration, assembly, & checkout 



7.5% x £(1) thru (14) 

16. 

Center support 

1 


15% 

x D(l) thru (15) 

18. 

Total spacecraft acquisition 

TACQC4 

N c 

x N /N . x 


L(l) thru (17) 


B. Research and Development 


1 . 

Prime power 

RDC 

2. 

Secondary power (batteries) 


3. 

Power conditioning 


4. 

Power distribution 


5. 

Antenna 


6. 

Transmitter 


7. 

Multiplexer 


8 

Receiver 


9. 

Structure 


10. 

Thermal control 


11. 

Stationkeeping 


12. 

Attitude control 


13. 

Telemetry and command 


14. 

Manned provision (life support subsystem) » 


(2) 


Curve 


Curve 

Thruput 

Curve 

i 

Curve 

Thruput 

Thruput 
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Table 2-1. Cost Model (continued) 

B. Research and Development (continued) 


15. 

Prototype 

RDC4(15) 

rACQC4(2) thru (14) 
+ 15% x ACQC3(1) 

16. 

Design, integration, and management 

RDC4(16) 

50% x £(1) thru (14) 

17. 

Center support and development 

RDC4(17) 

5% x [£(1) thru (14) +(16)1 

18. 

Ground support equipment 

RDC4(18) 50%x ^ACQC4(1) thru(14) 

19. 

Total spacecraft R/D 

TRDC4 

Z(l) thru (18) 


IV. LAUNCH VEHICLE 

A. Acquisition Cost 
Unit Acquisition 
Total Acquisition 

B. . Research and Development 
R/D Cost 

Total R/D Cost 
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N — total system lifetime in years 

L 

N — number of identical satellite systems 

s 

N — total satellite lifetime in years 

XjS 

N — number of satellites per launch 

SL 

Also N through N and M. through M 10 are variable and are included in the input 
, . , 1 16 1 12 
lists. 


The thruput cost values indicated for the launch vehicle are calculated external to 
the program and include the following items: 

Non-recurring : 

1) Launch vehicle development, 

2) Fairings and adapters, and 

3) Mission integration. 

Unit- recur ring : 

1) Launch vehicle, 

2) Fairings and adapters, 

3) Launch services, and 

4) Spares . 

2.5 ORBIT AND COVERAGE MODEL 

The model computes the slant range and elevation angle to the satellite and the off- 
axis location loss for a nominal ground station. These parameters are used in the 
propagation and noise models . 

To yield better results, the approach during the study has been to make use of a 
separate model, external to the program, prior to running a particular case. This 
program (see Figure 2-7) accepts the position of the satellite, the coordinates of 
the antenna beam center, the azimuth angle of the antenna, the antenna orthogonal 
beamwidths, and the relative angular positions of the desired radiation levels (e.g. , 
half-power, first sidelobe). Up to ten levels may be handled by the program. The 
program then computes the contours on the surface of the earth and plots on the 
SC-4200 plotter together with a mercator projection of the earth as in Figure 2-8. 
These contours are useful in determining the effectiveness of the coverage and in 
computing the flux density levels due to transmission. The model in the synthesis 
program accepts the same parameters to compute the elevation angle, major and 
minor axes , and location losses for each beam. 
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• SATELLITE POSITION 



• BEAM CENTER COORDINATES 

• BEAM AZIMUTH 


# ANTENNA PATTERN 


Figure 2-7* Orbit and coverage model. 



-too. Cl 35-0 -100.0 


0.0 3.0 X 7.0 


1070X770 


Figure 2-8. Continental United States coverage. 
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This section presents the various attenuation and noise factors which affect the pro- 
pagation between the ground uplink station, the satellite, and the ground receiving 
station and which are incorporated into the program. Associated curves are in the 
Appendix section of Volume 3. 

The basic model is shown in Figure 2-9. Given the elevation angles and the slant 
ranges for each satellite antenna beam, the attenuation due to clouds, rain, water 
vapor and oxygen, and the ionosphere as well as circuit losses are determined. 
Additionally, the noise contributions due to these elements are calculated. Table 
2-2 itemizes the terms and equations used in this model. The effective noise seen 
by the receiver is given by 


G 


T -f- T ■+■ T 
ANT T R 


where 


ANT 


= L 


TL 


(T, 


COS 


+ T + T 
I A 


+ T 


T ) 
RN ; 


G_ 


+ T G 
M M 


T G 
E E 


G , G G^ are relative antenna gains over the three regions — sky, man-made 
noise, and earth. 


The computation of the intermodulation noise for multicarrier operations is also 
made in this model. The relationship between carrier-to-intermodulation noise ( jyj) 
and transponder backoff for various numbers of carriers is shown in Figure 2-10. 
For a given c , or determined from the , the required transponder backoff is 

computed as a^function of the number of carriers . 

2.7 POWER AND SPACECRAFT SIZING SUBMODEL 


The power subsystem model is depicted in Figure 2-11. The model incorporates 
the various information types — audio, video, facsimile, and digital data as shown 
in the figure. Once the transmitter powers required to fulfill the links are deter- 
mined, the power requirements are traced back to the prime power source con- 
sidering the efficiency of each element in the system. 

The model will allow any number of channels of each type to operate at any duty 
cycle, and will size the power source and batteries for operation during both the 
sunlit and eclipse portions of the orbit. Each transmitter and receiver will operate 
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Figure 2-9. Synthesis program noise and propagation. 



BACKOFF FROM SINGLE CARRIER 
SATURATION -dB 


Figure 2-10. Transmitter backoff for multi -carrier operation. 
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Figure 2-11. Power subsystem model 
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at its full rate power or in a standby condition (nominally 10 percent full power) 
according to its duty cycle. 

The maximum power input to each transmitter and receiver block is 


P 

a 


P 

v 


P 


f 



( 

( 

(■ 

(■ 



Vta. *7 lea 


tv 


7?tv V lev 


tf 


77 tf 77 lef 


td 


r?tdTUcd 




then for maximum rating 


P =P.+P + P +P„ + P 

peu 1 a v f d , 


where P, = P ^ + P +P +p 

1 att stk mp ttc , 


P P. + P+ P+ P+ P 

_ peu _ 1 a v f d 


pc '/peu 


/peu 


now P 


P, + P* + P’ + P' + P' 
1 a v f d 

^ pc ^7 peu 


and P 


pti 


P+P + P +P„ + P, 
1 a v f d 

^pc^pcu 
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The total energy output from the batteries in watt-hours is 


H*P, + WH + WH + WH, + WH . 
WH . = 1 a v f d 

out 


where 


^pc^pcu 

H is the total number of hours per day of solar eclipse (1. 2 hours for a 
geostationary orbit) 


WH, 


= E (M-hj) [, 


I” BLOCK "1 ^ 


[ BLOCK "1 

['ON' powe rJ j + ^ 

( h j) 

[standby power] 


k = a, v, f, d 


h. is the number of hours during eclipse for standby of the j 
•* transmitter-receiver block 

The total energy supplied to the batteries is 
WH 


.th 


WH. = n 
in n i 


out 


HP. + WH + WH + WH. + WH 
1 a v f d 


^b^pc^pcu 


and the required power is 


WH. H*P, + WH + WH + WH, + WH^ 
_ in_ _ 1 a v f d 

pt 2 = 24-H ' (24-H)>?b>7pc>Jp c u 


The batteries are sized assuming a 50 percent depth-of-discharge requiring a 
capacity twice the output energy, WH out * 

Now the output of the power transfer is 

P + P +P + P+P. H*P +WH +WH +WH+WH, 
1 a v f d 1 a v f d 

pt pt^ ptg ^ pc^pcu (24 - H)7? b >7p C T?p CU 
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and finally the required power source output is 


P ^ - 


P s ^pt ^pt^pc 


~ — fp 

c”pcu J_ 1 


+ p +p +p + p + 
a v f d 


H*P +WH +WH + WH, + WH 
1 a v f d 


(24 - H)*7 b 


Because of the interdependence among the various satellite subsystems, an explicit 
solution of attitude control and stationkeeping parameters is precluded. The attitude 
control weight, volume, and power are a function of the solar array area and the 
stationkeeping characteristics are determined by the total of all of the remaining 
subsystem weights. Due to these relationships, an implicit solution is dictated as 
shown in Figure 2-12. The parametric values for each subsystem which is inde- 
pendent are determined and remain constant throughout the process. A nominal 
estimate is made for the parameters of the other subsystems — attitude control, 
stationkeeping, thermal control, and structure. The power subsystem is sized for 
these initial conditions. The structure is then sized as a function of the contained 
equipment; the thermal control requirements are determined for the required power 
dissipation. The attitude control subsystem parameters are computed for the solar 
array area and the stationkeeping subsystem is sized as a function of the weight of 
all other subsystems . The process is continued iteratively using the present state 
of values as the nominal until the difference in successive prime power require- 
ments is sufficiently small. 

2 . 8 REDUNDANCY MODE L 

The program computes factors which scale the affected weight, volume and cost 
terms to reflect the redundant elements required for additional spacecraft life. 

Let 

t = life of the satellite 
£ = life of the equipment in question 
P = a scale factor such that 0 < p < 1 

If d = — j- , then 

the scaling factor S is given by S = 1 + P • d. 


The same scale factor is used for weight, volume, and cost terms. 
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Figure 2-12. Iterative process for sizing satellite subsystems. 
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The list of equipment for which redundancy is provided includes the following: 

a) downlink transmitter 

b) spacecraft receiver 

c) power transfer unit 

d) power conditioning unit 

e) power control unit 

f) L/C filter 

Each requires that the lifetime l and the scaling factor P be input. 

2.9 COMMUNICATION LINKS 


The communication links are modeled using the relations shown in Figure 2-13 
where 


G 

is 

Antenna gain 

P 

is 

Transmitter power 

L 

is 

System attenuation 

N 

is 

System noise 

_C 

N 

is 

Carrier-to-noise ratio 

A 

is 

Satellite transponder gain 


Indicated subscripts are defined as follows 
U indicates 


D 

G 

S 

1 

2 


indicates 


uplink (at frequency f ) 

downlink (at frequency f ) 

ground based 

Satellite 

Station 1 

Station 2 


The receiver carrier-to-noise ratio of the signal transmitted by Station 1, trans- 
posed by the satellite, and received by Station 2 for a constant gain transponder is 


C_ 

N 


PG LG AG LG 
G UG U US DS D DG 

X 1 Li 

N~ + N„, + N A G~ L i” G„„ 

G IM S DS D DG 

Li Li 
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Figure 2-13. Communications link analysis 








GDC-DCF70-003 

where N is the intermodulation noise due to multicarrier operation through one 
IM 

transponder. 


For the constant power transponder, the equation relating satellite transmitter 
power and carrier-to-noise ratio is 


_C 

N 


(P 0 - A N ) L G^ _ 
v S S' DS D DG„ 


„ + N + N A G^ 0 L G„„ 

2 G IM S DS D DG 

Z Z 


In both equations, the total noise is composed of three terms: receiving system 
noise, intermodulation noise, and noise on the downlink due to the uplink. 


For the reverse link — station 2 to station 1 — the numerical subscripts are 
reversed . 


2.10. SELECTION OF DEPENDENT AND INDEPENDENT DESIGN PARAMETER 
VECTORS 


In this application, the choice of parameters for the dependent and independent 
vector varies with the case under consideration. Essentially, there are nine cases 
which the program will consider. The distinct cases are determined by the choice 
of transmit and receive capability of the Class I and/or Class II facilities and 
whether or not any direct receiving facilities are to be used. Figure 2-14 depicts 
the nine options which can be considered. Modifications to increase this number 
can be done with minimum effort. 

The choice for the components of Y, the vector of dependent parameters, is given as: 

y(l) — total satellite weight (lb) 

y(2) — total satellite volume (cu ft) 

y(3) — volume of the satellite equipment module (cu ft) 

y(4) — satellite prime power supply (watts) 

y(5), y(6), . . .y(28) — spacecraft transmitter power, power out/channel (watts) 
y(29), y(30), . . .y(52) — Class I transmitter power, power out/channel (watts) 
y(53), y(54), . . .y(76) — Class II transmitter power, power out/channel (watts) 
y(77) — direct ground antenna diameter or ground antenna diameter for 

selected dependent facility (ft) 
y(78) — direct ground antenna diameter (ft) 

For the components of the independent design parameter vector, the choice is: 

x(l) — Class I or Class II ground facility antenna diameter (ft) 
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Figure 2-14. Synthesis program ground system options 
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x(2) — Class I, Class II, or direct ground antenna diameter (ft) or Class I 

receiver noise figure 

x(3) — Class I, Class II, or direct receiver noise figure 

x(4) — Class II or direct receiver noise figure 

For the nine cases, the vectors are illustrated in Figures 2-15 and 2-16. The 
notation used in the charts is given by the following: 

W t = total satellite weight 

= total satellite volume 

Vem = volume of satellite equipment module 

P = prime power supply (satellite) 
ps 

Ptr = transmitter power (satellite) 

s 

Pg^ = transmitter power (Class I) 

Pg = transmitter power (Class II) 

A 

= antenna diameter (Class I) 

D = antenna diameter (Class II) 

A 

Dg = antenna diameter (direct) 

Nf^ = receiver noise figure (Class I) 

Nf = receiver noise figure (Class II) 

A 

Nf = receiver noise figure (direct) 

O 

In cases six through nine the choice of the independent antenna diameter is made on 
the basis of which class requires the largest value for C/N. The remaining antenna 
diameter reverts to being a dependent parameter. 


2.11 SOLUTION APPROACH 


As previously mentioned the link equations are based on the required ^ and this, 
in turn, is used to size the spacecraft. The key equations used in the power model 
are listed below. A more detailed discussion is found in Volume H which elaborates 
on the link model used in the program. 
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Figure 2-15. Choice of parameters for dependent vector, 
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CASE 

OPTION 

x(l) 

x (2) 

x(3) 

x(4) 

1 

D 

D 

Nf 



1 

3 

3 


2 

D 

D 

Nf 



1 

2 

2 


3 

D 

D 

Nf 

Nf 


1 

2 

2 

3 

4 

D 

1 

n, i 

— 

— 

5 

D 

Nf 

Nf 



1 

1 

3 


6 

D or D 

Nf 

Nf 



1 2 

1 

2 


7 

D or D 

Nf 

Nf 

Nf 


1 2 

1 

2 

3 

8 

D or D 

Nf 

Nf 



1 2 

1 

2 


9 

D or D 

Nf 

Nf 

Nf 


1 2 

1 

2 

3 

_ 


Figure 2-16. Choice of parameters for independent 
vector, X. 
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where R 


< G d>n < T j?d> nl 

<°d>k < T id>kj 


and the subscripts indicated are described as follows: 
i = data type = l, 2, 3, 4 

j = beam number = l, 2, 3 6 

k = fixed, representing independent class = 1 or 2 
n = fixed, representing dependent class = 1 or 2 



ratio of intermodulation noise/system noise 

ratio of satellite noise/receiver noise 

total noise at the ground receiver 

ground antenna downlink gain for class k 

spacecraft antenna downlink gain 

total loss down 

total noise at the spacecraft 

ground antenna uplink gain 

spacecraft antenna uplink gain 

total loss up 


Total noise on the ground, Ng, is a function of the effective temperature. The 
effective temperature at each of the ground receivers is a function of the relative 
antenna discriminations due to sky, environment, and earth. These discriminations 
in turn are determined by the downlink gain at each of the receivers. The gain is 
computed directly as a function of diameter and frequency, except when the diameter 
is unknown and is given as one of the dependent parameters. In this case, an itera- 
tion scheme is used since the equations are of an implicit nature. The iteration 
scheme is initiated using the value of the gain for the independent class of stations. 
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The steps describing the approach for each of the case options are given here. 


CASE OPTION I 

x - (D 2 . d 3 , n, 3 ) 

1. compute satellite antenna characteristics (from frequency and area of 
coverage), 

2a. compute Class I uplink gain, G , 

2b. compute direct antenna downlink gain, G^, 

3. compute antenna discrimination for Class I and direct, 

4. compute Class I uplink transmitter power, Pg^ and satellite transmitter 
power, Ptr, 

5. cost ground facility, 

6. size spacecraft and compute associated parameters (weight, volume, 
power, and cost), 

7. cost launch vehicle, 

8. compute total cost. 


CASE OPTION II 

* = < D i’ D 2’ Nf 2> 

1. I - 1, 

2a. I - 2a, 

2b. compute Class II antenna downlink gain, G , 

d u 

3. compute antenna discrimination for Class I and Class II, 

4. 1-4, 

5. I - 5 to I - 8. 


CASE OPTION III 

X = <Dj. D 2 , Nf 2 , Nf 3 ) 
i. ii - i to n - 4, 
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2. iterate for direct antenna diameter and gain, 

3. compute direct antenna relative discriminations, 

4. I - 5 to I - 8. 

CASE OPTION IV 

X = (D^ Nf x ) 

1. I - 1 and I - 2a, 

2. compute Class I antenna downlink gain, , 

3. compute Class I antenna discrimination, 

4. 1-4, 

5. I - 5 to I - 8. 

CASE OPTION V 

X = (D^ Nf^ Nf 3 ) 

1. I - 1, 

2. compute Class I antenna downlink gain, G „, 

dl 

3. compute Class I antenna discriminations, 

4. compute satellite transmitter power to be used in iteration scheme, 

Ptr (I, J) where I = JCLG and J = ffiMG, 

5. iterate for direct antenna diameter, D , 

6. I - 2a, 

7. compute direct antenna gain, G, . 

d3 

8. compute antenna discrimination for the direct antenna, 

9. I - 4 to I - 8. 

CASE OPTION VI 

X = (D , Nf^, Nf^) where j = 1 or 2 

1. I - 1, 

2. determine independent antenna diameter, G^ , and dependent antenna 

diameter, G. . 

dD 

3. compute downlink antenna gain for the independent antenna, 

4. compute antenna discriminations for the independent antenna, 
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5. V - 4, 

6. iterate for dependent antenna diameter, 

7. I - 2a, 

8. compute downlink gain for the dependent antenna, 

9. compute the antenna discriminations for the dependent antenna, 

10. I - 4 to I - 8. 

CASE OPTION VII 

X = (D., Nf , Nf , Nf ) where j = 1 or 2 

j 1 2 3 

1. VI - 1 to VI - 9, 

2. 1-4, 

3. Ill - 2 to III - 4. 

CASE OPTION VP! 

X = (D., Nf , Nf ) where j = 1 or 2 

J 1 ^ 

1. VI -1 to VI -9, 

2. compute Class II uplink antenna gain, G , 

3. 1-4, 

4. compute Class n uplink transmitter power, Pg Q , 

5. I - 5 to I - 8. 

CASE OPTION IX 

X = (D. , Nf , Nf , Nf ) where j = 1 or 2 

3 1 2 3 

1. vm-itovra-5, 

2. Ill - 2 to III - 4. 
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2.12 OUTPUT 

The ITS program provides two forms of output, one a brief summary printed at 
each iteration, and the second, a comprehensive printout at the conclusion of each 
case. 


A sample of the first type is illustrated in Figure 2-17. This summary is provided 
for each iteration and provides the program user with indication of the nature of the 
system at any particular design point during the convergence process and the rate 
of convergence. The output provides the iteration number, the values of the key 
dependent parameters, the values of the independent design parameters, the com- 
puted perturbation, and the gradient vector. It also provides the number of bound- 
ary violations (NBV), the index of the one which violates its constraint by the 
largest degree (KBV), the number of successive iterations for which the convergence 
criteria is met (ICV) , the number of components of X meeting the convergence 
criteria (JCV), the computed step size vector for the next iteration, and total 
system cost. 


The variables listed are defined as follows: 


WSAT 

VSAT 

VEM 

PPS 

D(l) or 

D(2) 

PG1 

PG2 

PTR 


total satellite weight (lb) 
total satellite volumes (cu ft) 

volume of the spacecraft equipment module (cu ft) 
prime power supply (watts) 

dependent Class I or Class II ground antenna diameter. 

4x6 array of transmitter power out at Class I facilities 
(watts) 

4x6 array of transmitter power out at Class II facilities 
(watts) 

4x6 array of transmitter power out at the spacecraft 
(watts) 


D(3), the direct receiving system antenna diameter is also provided when used. 


The formal output at the conclusion of the run is intended to be a comprehensive 
easy -to-read listing of all pertinent parametric values. 


Other than the output and computed values relevant to the user requirements and 
beam characteristics described in the previous sections, the output routine pro- 
vides a systematic breakdown of other elements — primarily cost, weight and 
volume. 
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Figure 2-17. Iterative output summary. 
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A total system cost summary is provided on an annual cost basis for each of the 
major subsystems: the Class I, Class II, and direct ground systems, the satellite, 
the launch vehicle, and the telemetry and command. 

A data service summary is provided for each of the four data categories on a per 
beam basis. This provides such information as base bandwidth, modulation, modu- 
lation index and total of bandwidth. (Figure 2-18) 

The ground receiver characteristics are listed for each class of stations and pro- 
vide a breakdown for each beam requiring that particular class . The information 
provides the number of channels transmitted and received, EIRP, G/T, S/N, and 
the total number of facilities for that beam. 

The satellite antenna and transponder characteristics are given for each beam. 
These include the number of channels transmitted and received, G/T, and EIRP. 
The transponder characteristics list the antenna dimension and gains, beam center, 
beamwidth, and pattern axes for the uplink and downlink. (Figure 2-19) 

A summary of the loss and noise terms are provided for the uplink and downlink 
beams. The noise terms include sky, earth, and environment, and receiver noise 
temperatures. The losses provided include free-space, atmospheric, rain, loca- 
tion, polarization, circuit, and pointing losses for the receiver and transmitter. 

On the downlink those applicable are given for each class of station. (Figure 2-20) 

Finally a comprehensive summary is provided for each of the major subsystems: 
ground, spacecraft, and launch vehicle. 

The ground facility summary is provided for each class of station and for each 
beam. The cost breakdown is given for each cost item comprising the ground 
system: facility, terminal equipment, transmitters, antenna and multiplexer, 
receiver, standby power, test equipment, personnel, and design, integration and 
management. The acquisition, installation, annual operations and maintenance, 
and total cost are listed for each. Also shown are the number of channels per data 
category, the required transmitter power, the antenna size and gain, and receiver 
noise figure. The total cost/facility is shown together with the total number of 
facilities for that particular beam. (Figures 2-21 and 2-22) 

The satellite summary provides the acquisition cost, R/D cost (prorated), total 
cost, weight, and volume of the various spacecraft subsystems: prime power, 
secondary power, power conditioning, power distribution, antenna, transmitter, 
multiplexer, receiver, structure, thermal control, attitude control, stationkeeping, 
and telemetry/command. Descriptive parameters such as transmitter power and 
receiver noise figures are provided where pertinent. Also included is the cost 
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Figure 2-18. Final output — user requirements /major subsystems costs. 
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satellite antenna characteristics 
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Figure 2-20. Final output — loss and noise summary 
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Figure 2-22. Final output — ground facility cost summary (2) 
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breakdown for integration, assembly and checkout, design, integration, and 
management, center support and ground support equipment. The total/satellite is 
given for the cost, weight, and volume. (Figure 2-23) 

The launch vehicle summary provides the acquisition, R/D (prorated), and total 
cost for the launch vehicle. Also shown is the weight and volume limits. 

A further breakdown for each beam is listed for the satellite communication sub- 
system. The summary provides cost, power, efficiency, lifetime for the trans- 
mitter, multiplexer, receiver, and antenna and the number of multiple transmitter 
blocks employed in the satellite. 
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CCSTS (Mill I CHS) HEIGHT VOLUME 

SATELLITE SUBSYSTEMS ACQUISITION *♦(> TOTAL lies.) (CU F T I TYPE DESCRIPTIVE PARAMETERS 

POWER SIBSYSTEN 

PRIME POWER 6. 15 3.16 11.33 *327.1 722.62 SOLAR POWER = 29. *2* KW 
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Figure 2-23. Final output — satellite subsystems and launch vehicle summary 
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APPENDIX A 

DERIVATION OF STEEPEST DESCENT ITERATIVE ALGORITHM 

Let X represent a vector of n independent system design parameters, 
y(X) a vector of m constrained system parameters, and C(X) the total system 
cost. It is desired to determine the minimum value of C subject to the 
constraints Imposed on y. This solution is to be obtained by iteratively 
perturbing a nominal design parameter vector X.^ until no father reduction 
in C is possible and all constraints are satisfied. At each iteration only 
those constraints which have been violated need be considered. 

Now let y be a k-vector whose elements are those parameters whose boundaries 
have been violated (a subset of the original y vector). For these parameters 
the constraint 

y(x) - y B0UMD = 0 (1) 

will be imposed. At a given iteration 
Y i = y ( X i) " y B0UND ^ ° 

Let g and <S> represent the cost and constraint gradients, respectively. That 
is, 

dC a 

« = dX ( X i) and g (X i ) (3) 

To first order terms 

ac(ax) = g T • AX (4) 

and AY = <5 ^ • AX ( 5 ) 

It is desired to determine the vector AX which minimizes 

J = AC(AX) = g T * AX ( 6 ) 

while reducing AY by seme specified amount 

AY - - K y • Y ( 7 ) 
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The magnitude of the perturbation is constrained by 

s 2 = AX T • uX (8) 

Applying Lagrange multiplier theory one defines J' to be 

J' = g T - AX + \ • (AX T • AX - S 2 ) + T1 T * ($ T • AX + K • Y) (9) 

V 

where \ is an arbitrary constant and T) is an arbitrary k-vector. From 
Equations 7 and 8 it is clear that J and J' are equal, so one can proceed 
by minimizing J ' . This minimum occurs where 


^ = g T + 2 X • AX T + if $ T = 0 


( 10 ) 


Solving Equation (10) for AX gives 


AX = - ^ • (g + $ * *H) (11) 

For convenience define 

K c = - and y = - ^ * T\ (12) 

Then 

AX=-K c *g-$-v (13) 

Substituting Equation (13) into Equations (5) and (j) gives 

AY = - K c <P T g - y = - K * Y (lb) 

Solving for y gives 

Y = - K c (<t T <b) • $ T • g + K y ($ T <p) • Y (15) 

Substituting Equation ( 15 ) into Equation ( 13 ) gives 


AX = - K [ I - $ (<3> T <p) $ T ] -g _K • $ ($ T $) • Y (16) 

y 

which is the perturbation equation given in Section 2.2. 
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The two terms of Equation (l6) represent perturbations parallel and perpendi- 
cular to the constraint boundaries, respectively. This can be shown by taking 
the vector inner product (vector projection) of the first term and the 
constraint gradient 




= - K c [ $ T - (<1? T $) (<P T $) $ T ] g 

= - K c [<P T - $ T ) g = 0 (IT) 


The value of K c can be determined by substituting Equation (l6) into 
Equation (8) and reducing to 


S 2 = K 2 Y T (A) Y + K 2 • H T H 
y c 


where 


H 


= [ I - <% £<*>) 1 


g 


(18) 

( 19 ) 


Then 


The step-size S is controlled to insure rapid convergence to the solution. 
When the nominal value of X is "far" frcm the solution successive perturba- 
tions will tend to be in the same direction. Under this condition it is 
desirable to increase S for the following iteration. When the nominal X 
is "close" to the solution such that a given perturbation overshoots the 
solution, successive perturbations will tend to be in opposite directions. 
Under this condition it is desirable to reduce S to allow the minimum cost 
solution to be determined to a finer resolution. The logic for implementing 
this step-size control is given in Section 2.2. 
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APPENDIX B 

GENERAL PURPOSE SUBROUTINES 

The ITS computer program requires the following general purpose system routines . 
Listings for each are provided. For the IBM 360 version, some may be replaced 
by the appropriate subroutines from the IBM library. 

General Purpose Routines 

1. Subroutine DSHIFT — to move an N by M matrix A into B 

2. Subroutine DZERO — to clear an N by M array 

3. Subroutine MTMPY — to form the product of two matrices 

4 . Subroutine MADD — to form the sum of two matrices 

Functions 

1. DVMAG — to compute the magnitude of a vector 

2. VDOT — to form the vector dot product of two vectors 

1. Subroutine DSHIFT 

Purpose: to move an N by M matrix A into B 

Usage: CALL DSHIFT (N, M, A, B) 

where A and B are assumed to be matrices dimensioned N by M 

2. Subroutine DZERO 

Purpose: to clear an N by M array 

Usage: CALL ZERO (A, N, M) 

where A is the array to clear and N and M are the dimensions of 
A 

3. Subroutine MTMPY 

Purpose: to form the product of two matrices C = A * B 
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Usage: CALL MTMPY (A, B, C, N, M, L) 

where A is an N by M matrix 

B is an M by L matrix, and 
C is the resultant N by L matrix 

Note: A and B may be the same matrix, but C must be a 
separate matrix. 

4. Subroutine MADD 

Purpose: to form the sum of two matrices R = F A + F A 

11 z z 

MADD (N, M, A x , F^ A^ F 2> R) 

A and A are N by M matrices, 

1 z 

F and F are scalar factors of A and A , and 
1 Z 1 z 

R is the resultant N by M matrix. 


Purpose: to compute the magnitude of an N dimensional vector V 
Usage: VM = VMAG (V, N) 

where V is an N-dimensional vector and VM is the resultant 
magnitude 

6. Function VDOT 

Purpose: to compute the scalar product of two vectors A and B 
S = A • B 

Usage: S = VDOT (A, B, N) 

where A and B are N dimensional vectors and S is the resultant 
scalar product 


Usage: CALL 

where 


5 . Function DVMAG 
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The following library functions 

FUNCTION 
SIN (X) 

ASIN (X) 

COS (X) 

AC OSD (X) 

TAN (X) 

A TAN (X) 

SQRT (X) 

A LOG 10 (X) 

EXP (X) 


are referenced: 

DEFINITION 

Sine x (radians) 

Arcsine x (radians) 

Cosine x (radians) 

Arccosine x (degrees) 

Tangent x (radians) 

Arctangent x (radians) 

Square root of X 

Log to the base 10 of X 
til 

e to the X power 
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SUBROUTINE DSHIFT ( I » J » A » B ) 
IMPLICIT REAL*B (A-H/O-Z) 
DIMENSION A(I,J),B(I»J) 

DO 1 K — 1 * I 
DO 1 L=1»0 
1 A(K»L) = B ( K ► L ) 

RE TURN 
EnU 


SUBROUTINE DZERO(A,I»J) 
IMPLICIT RtAL*B (A-HrO-Z) 
DIMENSION A ( 1 ) 

K — 1 * J 
DO 1 L=1»K 
1 A ( L ) = 0.0 
RETURN 
tNU 



SUbROuT I NE MTMPY ( A , B * C r M , N » L ) 

SUBROUTINE TO FORM THE PRODUCT OF TWO MATRICES 
A Is AN M BY N MATRIX 
B lb AN N BY L MATRIX 

C IS the resultant M BY l matrix 
implicit real*8 <a-h»o-z) 

DIMENSION A ( M » N ) > d(N,L), C(M,L) 

Do 20 1 = 1 »M 
DO 20 J— 1 » L 
C ( I 9 0 ) — 0.0 

DC 20 K-I » N 

20 C(IfJ) = CU»J) + A(1»K)*B(K#J) 

RETUR .. 

End 


subroutine madd (n » m» ai , fi * as»F2 »r ) 

SUBROUTINE TO ADD TWO N BY M MATRICES 
Al AND A2 ARE N BY M MATRICES 
R IS THE RESULTANT MATRICES 
Fl AND F2 ARE SCALAR FACTORS OF Al AND A2 

implicit real*** (a-h*o-z) 
dimension akd* A2(d» r(i> 

Ml - m*N 
DO 10o 1=1 » Ml 

10U R ( I ) = Fl * Al ( I ) + F2*A2(D 

K C turn 

END 
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Function dvmag<a,n> 
implicit real*8 (a-h»o-zj 
dimension a(1) 

Sum - o . 

DO 10 1=1»N 

10 SUM - SUM + A(I) * AtJ) 
OVMAG = DSURT(SUM) 

KETUKf- 

LNU 


FUNC I )ON V DO T ( A , o » N ) 
IMPLICIT RtAL*8 (A-H*0-Z) 
DIMENSION m(N), o(h) 

S = 0. 

DO 1 1 = 1 * N 

1 S = S + A(I ) * B ( I ) 

VDUT - S 

KLfUKfj 

END 
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APPENDIX C 

IBM 360 VERSION OVERLAY STRUCTURE 


Since the core storage requirements increase significantly in converting the pro- 
gram from the CDC 6400 version to the IBM 360 version, it is recommended that 
the overlay feature be used. The increased storage requirement, which could 
amount to 30% or more, is due to the shortened word length (32 bits/word vs. 60 
bits/word) and the need to revert to the use of double precision. 

Because of the fact that the major computations are contained within the do-loop 
system, the number of prospective subroutines to share the same levels is 
restricted; the philosophy of the overlay feature being to retain in core those mod- 
ules which are used consistently and flip in and out those which are used with less 
frequency. 

The resident level should then consist of the subroutines which are used for each 
iteration. This excludes the following list: 

Subroutine INPTX 
AOC 

CHANEL 

DIAM 

ATTEN 

PARBDS 

XBD 

These essentially are initialization routines and hence need not be maintained in 
core during the executions of the steepest-descent iterative process. For the CDC 
6400 version these represent approximately 13% of the allocated core. 

Subroutines OUTPT1 and OUTPT2 are used only once per case and need not be 
resident. These routines comprise approximately 41% of the core. 

The remaining ITS program subroutines and system subroutines make up the 
remainder; the driver routine requiring only 6% of the core requirement. There- 
fore to implement the overlay, three levels of routines may be considered: 
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Resident: 

Driver routine (6%) 


& level: 

Initialization routines (13%) 

INPTX 

AOC 

CHANEL 

DIAM 

ATTEN 

PARBDS 

XBD 

level: 

System routine (38%) 


EQUAT 

PXTR 

SPCRFT 

RCVRS 

OPTMUM 

DISCR 

ANTS 

STNKP 

BNDV 

DGDEP 

ATTCN 

STRUC 

CNVRGE 

GRDFAC 

MANPRV 

TLMCMD 

DIAGN 

GDIRCT 

MTPLXS 

THRML 

MAXBD 

LNCHV 

POWER 

XMTR 


TOTCST 


PLACE 


V level: 

Output routine (41%) 

OUTPT1, OUTPT2 

Added to these are the required operating system routines. 

To accomplish the overlay, only one level is considered — level oi. Three separate 
packages are interchanged during execution. 

a - INPTX 

1 

a — EQUAT, OPTMUM, BNDV / CNVRGE, MAXBD, PXTR, DISCR, 

2 DGDEP, POWER, RCVRS, STNKP, STRUC, TLMCMD, THRML, 
XMTR, PLACE, GRDFAC, GDIRCT, LNCHV, TOTCST, SPCRFT, 
ANTS, MANPRV, MTPLXS, ATTEN, VDOT, DVMAG 

a - OUTPT1, OUTPT2 

3 

All else is resident. 
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